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Abstract. We present a clear detection of CH2 in absorption towards the molecular cloud complexes Sagittarius B2 and W49 N 
using the ISO Long Wavelength Spectrometer. These observations represent the first detection of its low excitation rotational 
lines in the interstellar medium. Towards Sagittarius B2, we detect both ortho and para transitions allowing a determination of 
the total CH2 column density of 2V(CHj) = (7.5 ± 1.1) x 10 14 cnT 2 . We compare this with the related molecule, CH, to determine 
[CH/CH 2 ] = 2.7 + 0.5. Comparison with chemical models shows that the CH abundance along the line of sight is consistent with 
diffuse cloud conditions and that the high [CH/CH2] ratio can be explained by including the effect of grain-surface reactions. 
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1. Introduction 

Methylene (CH2) is thought to be a relatively abundant 
molecule in diffuse as well as in dense interstellar clouds, with 
similar abundances to CH (e.g. van Dishoeck & Black 1986; 
Lee et al. 1996). However, it has proved very difficult to detect 
observationally, mainly due to the inaccessibility of its rota- 
tional lines from the ground. 

Methylene was initially proposed to explain an unidentified 
ultraviolet (UV) band observed in comets (Herzberg 1942a,b) 
but later this band was shown to be associated with C3 (Douglas 
1951). CH2 remains undetected in comets, but recently, several 
of its electronic bands have been observed in the interstellar 
medium (ISM) by Lyu et al. (2001) using the Hubble Space 
Telescope. They tentatively detected CH2 absorption bands in 
the UV spectrum towards two stars, HD 154368 and £ Oph. 

Interstellar searches for rotational emission/absorption are 
hampered by the peculiar spectrum of CH2, caused by its light- 
ness and b-type selection rules (Michael et al. 2003). This re- 
sults in lines at widely varying wavelengths, all of which are ei- 
ther completely unobservable from the ground or are in difficult 
spectral regions that are at the edges of atmospheric windows 
or for which there are few telescopes with suitable instrumen- 
tation to observe them. 
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Hollis et al. (1995) have, so far, made the only unambigu- 
ous identification of CH2 in the ISM, confirming their earlier 
detection (Hollis et al. 1989). They clearly identified the 4o4- 
3 13 rotational transition with simultaneous measurements of 
multiple fine-structure features between 68 and 71 GHz. These 
were detected in emission towards Orion KL and W5 1 M, both 
dense "hot core" sources, which provide excitation to the 4o4 
level which is » 215 K above the ground state. The lines oc- 
cur in a rarely observed spectral region, close to the telluric 
55-65 GHz O2 bands. Accurate frequencies for this transition 
were measured by Lovas et al. (1983). 

Low-excitation rotational transitions occur in the far- 
infrared (FIR) region, accessible by the Infrared Space 
Observatory (ISO) satellite (and in the future by SOFIA and 
Herschel), and around 320 fim (940 GHz), a spectral region just 
being opened from the ground. Although limited laboratory 
data on the rotational spectrum of CH2 have existed for many 
years (Lovas et al. 1983; Sears et al. 1984), much more accu- 
rate frequencies are currently being measured at the Cologne 
Laboratory for Molecular Spectroscopy (Michael et al. 2003; 
Briinken et al. 2004) 

We have used these new data to search for low-lying rota- 
tional transitions of CH2 observed by the ISO Long Wavelength 
Spectrometer (LWS; Clegg et al. 1996). We clearly detect 
the strongest fine structure components of transitions from 
the lowest energy level of both ortho and para CH2 towards 
Sagittarius B2 (Sgr B2) and of ortho CH 2 towards W49 N. Both 
Sgr B2 and W49 are giant molecular cloud complexes emitting 
strong FIR continuum spectra. This makes them particularly 
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Fig. 1. Low-lying rotational levels of the CH2 molecule showing the 
wavelength of the strongest fine structure component for each transi- 
tion. 



good targets for detecting absorption lines from intervening in- 
terstellar matter. 

Sgr B2 was observed as part of a wide spectral survey us- 
ing the LWS Fabry-Perot (FP) mode, allowing us to search for 
transitions from all the other low-lying energy levels, as well 
as to compare the data with absorption from the chemically re- 
lated CH molecule. We use the strong CH lines as a template to 
fit and calculate column densities for CH2. In Sects. 4, 5 and 6 
we present the results towards Sgr B2 and in Sect. 7 detail the 
results for W49 N and several other sources where CH2 was not 
detected. We then compare the results with chemical models in 
Sect. 8. 



2. CH 2 rotational spectrum 

Figure 1 shows the low-lying rotational states of CH2. The en- 
ergy levels are denoted by N KaKc , where N is the rigid-body 
rotational quantum number. Each of these levels is split by 
electron spin-spin and spin-rotation interactions into three fine- 
structure levels, with the quantum number for total angular mo- 
mentum excluding nuclear spin, 7, equal to N + 1, N, N - 1 for 
N > (7=1 at N=0). Due to the presence of two protons with 
opposite nuclear spin, CH2 has both ortho (nuclear spin quan- 
tum number, 1=1) and para (7=0) forms. In the ortho levels, 
the non-zero nuclear spin interacts with the electrons to cause a 
further hyperfine splitting with total angular momentum quan- 
tum number, F=J + 1, 7, 7 - 1. Selection rules for rotational 
transitions are A7=0, ± 1 and AF=0, ± 1 . 

Frequencies for low-lying rotational transitions have been 
calculated from Laser Magnetic Resonance (LMR) laboratory 
spectra by Sears et al. (1984) with a quoted accuracy of 5 MHz 
(0.6 km s _1 ). In order to confirm the lines and to calculate ac- 
curate line strengths we have used new measurements of the 
sub-mm lines (Michael et al. 2003; Briinken et al. 2004) to 
determine more accurate frequency values. Table 1 shows the 
calculated wavelengths, Einstein coefficients and line strengths 



Table 1. Calculated wavelengths, Einstein coefficients (Ay) and 
line strengths (Sy) for the low-lying transitions of CH2 (ignoring 
hyperfine-structure in the ortho transitions). The estimated error corre- 
sponding to the last significant figure of the calculated wavelengths is 
given in parenthesis. The Einstein coefficients and line strengths cal- 
culated here are half those given by Sears et al. (1984). The detected 
lines towards Sgr B2 are shown in bold. 
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Wavelength 


Ay 
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1.14 
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153.814187 (6) 


0.00962 


5.15 


9, 
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152.992296 (6) 


0.00216 


1.14 




-2m 7=3-2 


154.178941 (6) 


0.00151 


1.14 


2 n - 


-2 02 7=3-3 


153.352914 (7) 


0.0125 


9.25 


3i 3 - 


-2 02 7=2-1 


93.5838 (1) 


0.0298 


3.60 


3i 3 - 


-2o2 7=2—2 


94.0309 (1) 


0.00559 


0.685 


3i 3 - 


-2 02 7=2-3 


93.7231 (2) 


0.000170 


0.0205 


3l3- 


-2 02 7=3-2 


93.7025 (2) 


0.0317 


5.37 


3l3- 


-2 02 7=3-3 


93.3967 (2) 


0.00392 


0.657 


3l3- 


-2 02 7=4-3 


93.6621 (2) 


0.0356 


7.75 


lio- 


-loi 7=0-1 


156.43947 (1) 


0.0133 


0.500 


lio- 


-1 , 7=1-0 


155.66632 (1) 


0.00462 


0.513 


lio- 


-loi 7=1-1 


157.56539 (1) 


0.00318 


0.366 


lio- 


-loi 7=1-2 


156.76825 (1) 


0.00548 


0.622 


lio- 


-loi 7=2-1 


157.12503 (2) 


0.00323 


0.615 


lio- 


-loi 7=2-2 


156.33233 (1) 


0.0101 


1.88 


2l2- 


-loi 7=1-0 


107.2947 (2) 


0.0134 


0.487 


2l2- 


-loi 7=1-1 


108.1935 (2) 


0.0103 


0.385 


2l2- 


-loi 7=1-2 


107.8170 (2) 


0.000731 


0.0270 


2l2- 


-loi 7=2-1 


107.8573 (2) 


0.0184 


1.13 


2l2- 


-loi 7=2-2 


107.4832 (2) 


0.00603 


0.367 


2l2- 


-loi 7=3-2 


107.7203 (2) 


0.0244 


2.10 



(averaging over the hyper-fine structure in ortho transitions). 
The accuracy of the calculated wavelengths is highest for those 
transitions actually measured in the lab (2n-2o2 (153 //m) and 
1 10— loi (156 yum); Briinken et al. 2004). The uncertainty in the 
calculated Einstein values and line strengths is dominated by 
that of the dipole moment. We have followed previous authors 
and used a value of 0.57 D (derived from an ab initio calcula- 
tion by Bunker & Langhoff 1983). The calculation strongly de- 
pends on the electronic and geometric structure of the molecule 
and inserting the most recent values for the geometry suggests 
that the uncertainty must be at least 0.02 D (leading to ~5% 
uncertainty in the Einstein coefficients). The line strength in 
Table 1 is defined in the same way as by Sears et al. (1984), al- 
lowing a direct comparison. This shows that the strengths cal- 
culated here are half those previously published. However, our 
results can reproduce the values for the 1 n-2 2 (317 //m) tran- 
sition given by Michael et al. (2003) and agree with the calcu- 
lations of Chandra (1984) (both of these authors used an inde- 
pendent method to calculate their line strengths). This indicates 
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Table 2. Coordinates of observed sources. 



Source 


RA (J2000) 


Dec. (J2000) 


SgrB2 


17 h 47 m 21.75 s 


-28°23'14.1" 


W49N 


19 h 10 m 13.55 s 


+09°06'14.7" 


Sgr A' 


17 h 45 m 39.97 s 


-29°00'23.6" 


NGC7023 


21 h 01 m 36.91 s 


+68°09'48.2" 



that the line strengths quoted in the paper by Sears et al. (1984) 
are likely to be in error by a factor of two (T. Sears, private 
communication). 

3. Observations and data reduction 

The Sgr B2 observations were carried out as part of a wide 
spectral survey using the ISO LWS Fabry-Perot (FP) mode, 
L03. The whole LWS spectral range, from 47 to 196 /J.m (6.38- 
1.53 THz), was covered using 36 separate observations with a 
spectral resolution of 30-40 km s _1 . The first results from this 
survey have recently been presented by Ceccarelli et al. (2002); 
Polehampton et al. (2002); Vastel et al. (2002); Polehampton 
et al. (2003). We also analysed additional L04/L03 mode ob- 
servations towards Sgr B2, W49 N, Sgr A* and NGC7023 and 
one grating L01 observation towards W49 N, downloaded from 
the ISO Data Archive 1 . Details of the ISO TDT numbers of all 
the observations used in this paper are shown in Table A. 1 . 

The LWS beam has an effective diameter of ~ 80" (Gry 
et al. 2003) and was centred on the coordinates shown in 
Table 2. The coordinates chosen for Sgr B2 were slightly offset 
from the main FIR peak at Sgr B2 (M) to ensure that Sgr B2 (N) 
was excluded from the beam. Three positions were observed 
towards NGC7023 - the illuminating star (coordinates given 
in Table 2), the NW PDR (-24.5", +39.2") and the SW PDR 
(-25.8", -63.1"). 

The L03 observations towards Sgr B2 were reduced using 
the ISO Offline pipeline (OLP) version 8 with the remaining 
observations processed with OLP version 10 (for FP observa- 
tions there is no significant difference between OLP 8 and 10). 
Further processing was then applied interactively using rou- 
tines that appeared as part of the LWS Interactive Analysis 
package version 10 (LIA10: Lim et al. 2002) and the ISO 
Spectral Analysis Package (ISAP: Sturm et al. 1998). The inter- 
active processing included the calculation of accurate dark cur- 
rents (including stray light), careful removal of the LWS grat- 
ing profile from the spectra and removal of glitches caused by 
cosmic ray impacts (see Polehampton et al. 2002, 2003). The 
wavelength scale of each observation was then corrected to the 
local standard of rest and data co-added for each line (small 
multiplicative factors were sometimes necessary to match each 
observation in flux before co-adding). The final uncertainty in 
wavelength is less than 0.004 //m (or 11 km s~') correspond- 
ing to the error in absolute wavelength calibration (see Gry 
et al. 2003). Finally, a polynomial baseline was divided into 
the data to obtain the line-to-continuum ratio. For the Sgr B2 

1 see http://www.iso.vilspa.esa.es/ida. 



data, a 3rd order polynomial was necessary to fit the contin- 
uum around and between the detected lines. For the L04 data 
towards W49 N where the continuum coverage is much lower, 
a 1st order baseline was used. The use of low order polynomi- 
als and careful masking of the detected lines for the fit ensured 
that no spurious features were introduced into the data by the 
division. This step in the reduction is very useful because it ef- 
fectively bypasses the large uncertainty in absolute flux level 
caused by multiplicative calibration steps (see Swinyard et al. 
1998). The remaining error is dominated by statistical noise in 
the data with a small additional uncertainty due to the dark cur- 
rent determination and continuum fit. 

In Sect. 6, we use one grating L01 observation towards 
W49 N to determine the column density of CH. This observa- 
tion was reduced using OLP 10, and corrected for saturation in 
the LWS detectors (necessary due to the strength of the W49 N 
thermal continuum) by T. Grundy at the UK ISO Data Centre. 
The dataset was further reduced using the LWS L01 post- 
processing pipeline which performs several additional correc- 
tions detailed in Lloyd et al. (2003). This allowed us to cal- 
culate an accurate value for the line-to-continuum ratio in the 
line. 

4. Sgr B2 results 

We searched the Sgr B2 spectrum for all the low-lying rota- 
tional lines of both ortho- and para-CH2 occurring within the 
survey range (47-196 /J.m). Figure 2 shows the data around the 
three lowest energy transitions (fits to the line shapes are also 
shown - these are described in Section 4.4). At the spectral res- 
olution of the LWS FP, the fine structure splitting of each ro- 
tational transition is resolved but the spacing of the hyperfine 
components in ortho-CPL. (< 10 km s~') is too small to be sep- 
arated. When analysing ortho transitions we have used wave- 
lengths averaged over the hyperfine structure weighted by A, ; -. 
The wavelengths used are shown in Table 1 . 

The survey also shows strong absorption due to CH (shown 
in the top panel of Fig. 2). We have used these lines as a basis 
for fitting the observed CH2 line shapes to fix the relative con- 
tribution from features along the line of sight. This was nec- 
essary because the signal-to-noise in the CH2 detections is not 
high enough to allow a full fit accounting for different line of 
sight components separately. The method assumes that there is 
a constant [CH/CH2] ratio applicable to all line of sight com- 
ponents. 

4.1. CH 

Two prominent absorption features due to CH are present in the 
spectrum due to its 2 ni/2 7=3/2-1/2 transition from the ground 
state to the first rotational level, ~96 K above ground (see also 
Cernicharo et al. 1999; Goicoechea et al. 2004). These are due 
to the A-doublet type splitting of each rotational level and oc- 
cur at wavelengths of 149.09 //m and 149.39 /im (Davidson 
et al. 2001). The two detected lines are shown in the top panel 
of Fig. 2 and are broad with absorption in the range -150 to 
+ 100 km s _1 . This is due to galactic spiral arm clouds be- 
tween the Sun and Galactic Centre (centred at velocities -100 
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to +30 km s -1 , e.g. see Greaves & Williams 1994) as well as to 
Sgr B2 itself. The peak absorption occurs near the velocity of 
Sgr B2 at ~ +65 km s _1 . These lines have previously been de- 
tected with the Kuiper Airborne Observatory (KAO) by Stacey 
et al. (1987) but the ISO data provide a significant improve- 
ment in the signal-to-noise ratio as well as complete coverage 
of the neighbouring continuum. No higher transitions of CH 
were detected in the survey above the noise. 

In order to determine the CH column density for each ve- 
locity component in the line of sight towards Sgr B2, a high 
resolution model of the line shape was constructed, convolved 
to the LWS resolution and adjusted to obtain a best fit. This 
was carried out in the same way as used for OH lines (see 
Polehampton et al. 2003, Polehampton et al. in preparation). 
H 1 21 cm absorption measurements (Garwood & Dickey 1989) 
were used to fix the velocity and line width of 10 line of sight 
components and then optical depths were adjusted in a multi- 
parameter fit to find the minimum in^ 2 . The optical depth for 
each component, t, was calculated from the line-to-continuum 
ratio using 

/ = / c exp(-r) (1) 

where I c is the intensity of the continuum. This method as- 
sumes that the same velocity components seen in the H i spec- 
trum are present in CH and CH2 (but the fit does not depend 
on the H 1 optical depths derived by Garwood & Dickey 1989). 
This is likely because the atomic material is seen to be asso- 
ciated with molecular gas at the same velocities (e.g. Vastel 
et al. 2002). Each fitted component probably still represents a 
mean over many narrower features such as those seen in CS ab- 
sorption with velocity widths ~1 km s _1 (Greaves & Williams 
1994). The final fit is shown in the top panel of Fig. 2. 

Column densities for each velocity component are shown in 
Table 3. These were calculated assuming a Dopplerline profile 
with Maxwellian velocity distribution (e.g. Spitzer 1978), 

= ^£ 10 " ^ (2 ) 

2VhTI AijA'jgi/gj 

where Nj is the column density in the lower level, to is the 
optical depth at line centre, Av is the line width in km s _1 , A, ; - 
is the Einstein coefficient for spontaneous emission, Ay is the 
wavelength in //m and gj is the statistical weight of state i. 

No higher transitions of CH have been observed in the spec- 
trum even at the velocity of Sgr B2 itself, showing that the 
level populations are characterised by a low rotation temper- 
ature (i.e. sub-thermal excitation) and the ground state popula- 
tion is a good measure of the total column density. Using a 2cr 
detection limit for the 2 ni /2 7=5/2-3/2 CH transition (1 16jum), 
we calculate r rot < 20 K. The fitted optical depths show that the 
absorption is generally optically thin, although in the Sgr B2 
component at 67 km s _1 it is marginally optically thick (t=2.4). 
The results of the fit are shown in Fig. 2 and Table 3, with er- 
rors determined by examining^ 2 as pairs of optical depths were 
varied about their best fit values. 

One way to check the opacity of Sgr B2 is to use the less 
abundant 13 CH isotopomer whose 2 ni/2 7=3/2-1/2 transition 
has components at 149.79 jim and 150.09 //m (Davidson et al. 



Table 3. Results of our fit of the CH lines towards Sgr B2 are given 
in column (A). The velocities and line widths are taken from the Hi 
data of Garwood & Dickey (1989). Column (B) gives the previous 
estimates of CH column densities from the KAO measurements of 
Stacey et al. (1987) and columns (C) and (D) give the results from 
the radio A-doublet lines from Genzel et al. (1979) and Andrew et al. 
(1978). 



LSR Velocity 


FWHM 


N(CH) (10 14 


cm 2 ) 




(kms- 1 ) 


(km s" 1 ) 


(A) 


(B) 


(C) 


(D) 


-107.6 


7 


0.9+0.7 


4.2 






-81.7 


28 


2.6+1.1 






2.5 


-51.9 


17 


0.3±0.15 








-44.0 


8 


1.5±0.6 


3.6 




1.3 


-24.4 


14 


1.9±0.4 






0.3 


+ 1.1" 


19 


2.0±0.3 


3.2 




1.1 


+ 15.7 


7 


3.1+0.8 






1.5 


+31.4 


21 


1.5+1.0 








+52.8/ + 66.7* 


11/16 


9.3±0.9 


5.2 


3.1 


9.4 C 


Total 




20.1 


15.2 




16.1 



" This velocity is corrected from an error in Garwood & Dickey (1989) 
where it is given as +1 1 km s -1 

* The Hi data resolve 2 components in Sgr B2. However, our fit re- 
quires only one of these (at 66.7 km s~') to reproduce the ISO spec- 
trum. 

c This is the combined value for all components > 50 km s -1 assuming 
the emission fills the beam, as described in the text by Andrew et al. 
(1978). 



2004). Langer & Penzias (1990) find a 12 C/ 13 C ratio towards 
Sgr B2 of 24 ± 1 using observations of 12 C ls O and 13 C 18 0. 
However, the 13 CH lines are not detected in the ISO spectrum 
to a limit of 1.3% of the continuum (2<x level). This leads to 
12 CH/ 13 CH > 36 at the velocity of Sgr B2. This is higher than 
the value derived from CO but this could be due to isotopic 
fractionation which is expected to increase the ratio from the 
true 12 C/ 13 C value (see Langer et al. 1984). In the following 
analysis, we use the column density for Sgr B2 without adjust- 
ment. 

The final CH column densities in Table 3 compare well 
with previous values derived from the KAO observations of 
same FIR lines by Stacey et al. (1987). Goicoechea et al. (2004) 
have also examined ISO data of CH and found absorption ex- 
tended across the whole surrounding region. They calculate to- 
tal column densities in the range (0.8-1. 8)xl0 15 ctrT 2 , peaking 
at the central Sgr B2 (M) and Sgr B2 (N) positions. Radio ob- 
servations of the CH A-doublet lines by Andrew et al. (1978) 
are in very good agreement with the individual velocity fea- 
tures that we have fitted - the small discrepancies may be due to 
their much larger beam size (8.2'), which also included the (N) 
position. The radio lines were also observed by Genzel et al. 
(1979), however, they only give the column density for Sgr B2 
itself at +65 km s _1 (although they also observed emission from 
the other velocity components). They derived the column den- 
sity from the 3264 MHz line giving a value three times lower 
than ours (however, their weaker data for the 3335 MHz line 
would give a higher column density). 
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4.2. Ortho-CH 2 

The lowest ortho level of CH2 occurs at the ground level, N-0. 
We clearly detect absorption from the 1 n-Ooo J -2-1 and 7=1- 
1 fine structure components at 127 /on (see Fig. 2). The third 
component, 7=1-0, is not detected above the noise in the spec- 
trum. The observed features have a similar shape to the CH 
lines, indicating the presence of absorption along the whole line 
of sight. 

The next highest energy level occurs ~67 K above ground. 
There are two transitions from this level that occur within 
the ISO wavelength range at 153 /im (2n-2o2) and 94 /mi 
(3i3-2o2), and one that occurs outside the range at 317 /mi 
(ln-2o2). The strongest fine-structure component of these is 
the 2n-2o2 7=3-3 line at 153.353 /mi. However, this wave- 
length is very close to the iV=2-l 7=3-2 transition of NH at 
153.348 /mi (separation ~10 km s _1 ) and there is a strong ab- 
sorption line detected (see Cernicharo et al. 2000). The identi- 
fication of this feature with NH is secure because we also de- 
tect its 7=2-1 and 7=1-0 lines with relative depths that agree 
extremely well with the predicted line strengths (measured val- 
ues are 1.0/0.55/0.27 compared to the predicted relative line 
strengths of 1.0/0.54/0.23 from the JPL line catalogue; Pickett 
et al. 1998). Furthermore, the next strongest CH 2 transition at 
153.814 /an is not detected. This shows that the major contri- 
bution to the 153.35 /mi line must be from NH rather than CH2. 

The next strongest transition is 3i3-2o2 and the 7=4-3 line 
(93.662 //m) should be almost as strong as the 153.35 /mi line 
discussed above. There appear to be some features at roughly 
the correct velocity at a level ~ 3cr above the statistical noise 
(see Fig. 2 but note that the error bars shown include systematic 
uncertainty). 

No higher level ortho transitions were detected above the 
noise in the data. 

4.3. Para-CH 2 

The lowest energy state for para-CH2 occurs at ~23 K above 
ground (see Fig. 1) and has two transitions at 156 /im (loi-lio) 
and 107 /im (2i2— 1 10)- The strongest fine-structure component 
of these is the 2i2-l 10 7=3-2 transition at 107.720 /mi. There is 
an absorption feature at approximately the correct wavelength 
for this line, although it appears to be wider than expected in its 
short wavelength wing. The observed absorption depth is also 
consistent with the data at the position of the next strongest 
component, 7=2-1, at 107.857 /im (see Fig. 2). Using this de- 
tection to predict the optical depth of the strongest loi-l 10 tran- 
sition at 156.332 /mi (7=2-2) gives a value consistent with the 
noise level in the spectrum. 

No higher para transitions were detected above the noise in 
the data. 



4.4. CH 2 column densities 

In order to fit the CH 2 lines, we assumed that the relative col- 
umn densities between each velocity feature were the same as 
for CH. This fixes the shape of the line and allows a single pa- 
rameter fit to determine an 'average' CH/CH2 column density 
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Fig. 2. Upper panel: Fit of the CH line towards Sgr B2 with the data 
binned at a quarter of the resolution element. Lower panels: In the 
lower three panels, a model based on the CH line shape has been fitted 
to the data around the three lowest energy transitions of CH2. The rel- 
ative contribution of each velocity component in these fits was fixed 
to the CH values in Table 3, except for the 94 /im (3 13 — 2 2) transi- 
tion in the bottom plot which only includes the velocity component 
due to Sgr B2 itself (see text). The fit was carried out to the strongest 
fine structure component in each case. Other (sometimes undetected) 
components are included in the model using their expected relative 
line strengths from Table 1 . The data are binned at half of the resolu- 
tion element and the errors reflect both statistical noise and systematic 
error due to the dark current subtraction and continuum fit. The extra 
absorption feature seen at -250 km s -1 relative to the 127 /mi (ln- 
Oqo) transition appears to be real but is not identified. 
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ratio for each CH 2 energy level. The fit was carried out to the 
strongest fine structure component in each case, with the rela- 
tive importance of the other (sometimes undetected) fine struc- 
ture transitions fixed using the line strengths from Table 1 . The 
results of these fits are shown in the lower 3 panels of Fig. 2 (in- 
cluding the undetected transitions, to show that their predicted 
absorption is consistent with non-detection in the data). 

This method worked well for the strongest transition (ln- 
Ooo) at 127 yum, where both the line shape and relative strength 
of the J -2-1 and 7=1-1 components show a good fit. A detec- 
tion limit for the 7=0-1 component is also consistent with the 
expected line strengths of 1.0/0.6/0.2 for 7=2-1, 7=1-1 and 
7=0-1 respectively. The fact that the line shape agrees well 
with the data indicates that the [CH/CH2] ratio does not vary 
significantly between the different line of sight components. 

The para transition at 107.7 //m (2i2— lot) does not show 
such good agreement in line shape with the model. As men- 
tioned in Sect. 4.3, the main line appears wider than expected, 
although the noise is high in this part of the spectrum. It is un- 
likely that this extra absorption is associated with the CH2 line 
as this would mean that the negative velocity features would 
have to be stronger than the component due to Sgr B2 itself 
and this is not observed in other lines. It could be due to overlap 
with an unidentified feature (although we have not been able to 
find a likely candidate) or an instrumental effect (e.g. the over- 
lap of mini-scans can sometimes cause problems for wide lines 
- see Polehampton et al. 2003). A final possibility is that the ex- 
tra absorption could have been introduced when dividing by the 
polynomial fit to the continuum. However, we minimised this 
effect by using a low order polynomial (3rd order) to carefully 
fit only the continuum around features that already appeared in 
the raw data. Any variation in the fitted baseline is on a larger 
scale than the excess absorption in the line. In fitting the line, 
we have only considered velocities above km s" 1 . The actual 
line shape is not well constrained by these data and our fitting 
method of fixing the relative contribution of each velocity to 
that of CH may overestimate the absorption due to the line of 
sight clouds. 

In order to fit the data around the transition at 94 //m (3 13- 
202), we modified the line template to include only the velocity 
component at +67 km s _1 from Table 3. This is because we 
only expect a significant population in the higher levels at the 
velocity of Sgr B2 itself. In Sgr B2 itself, excitation to the 2o2 
level (~67 K above ground) can be provided by FIR photons 
from the strong dust continuum, whereas in the line of sight 
clouds there is no strong radiation field and the densities are 
too low for collisional excitation to be important. The bottom 
panel in Fig. 2 shows the results of the fit. 

The small discrepancy in velocity between fit and data can 
be explained by the fact that the main feature is only 3cr above 
the noise. The uncertainty in the exact line wavelength is also 
relatively high for this transition as it has never been measured 
in the lab either in Cologne or by Sears et al. (1984). The dif- 
ference between the value quoted in Table 1 and that by Sears 
etal. (1984) is ~3kms-'. 

The final column densities for the lowest ortho and para lev- 
els of CH2 were determined directly from the fitted CH/CH2 ra- 
tio for each level, assuming that the distribution over fine struc- 



ture levels is determined by the line strengths given in Table 1 . 
The final best fit values summed over all velocities are, 

N(0 00 ) = (2.9 ± 0.3) x 10 14 cm" 2 
AT(loi) = (3.4 ± 0.9) x 10 14 cm" 2 

The column density in the 2o2 level at the velocity of Sgr B2 
is, JV(2 2) sg rB2 = (1.2 ± 0.5) xlO 14 cm" 2 . Treating this as an up- 
per limit allows us to make a comparison with the ortho ground 
state (using only the Sgr B2 component) and estimate the ex- 
citation temperature for Sgr B2, giving T mt <40^J 4 K. In or- 
der for absorption to be seen, this excitation temperature must 
be lower than the temperature of the dust producing the FIR 
background. In Sgr B2 the dust temperature has been found to 
be ~ 30 K (Goicoechea & Cernicharo 2001), and so our esti- 
mated excitation temperature is consistent within its errors. At 
this temperature we do not expect any higher levels to be sig- 
nificantly populated and do not detect any of the other higher 
transitions that occur within the spectral survey range (220— 1 1 1 
at 50.5 yum; 221— 1 10 at 50.8 //m and 3i2— 3o3 at 148.8.5 /mi). 

Assuming that only the first three energy levels in Sgr B2 
are populated, the ortho-to-para ratio in the +67 km sT l com- 
ponent is 1.6^'g. The equilibrium value of the ratio for CH2 
should be >3 for low temperatures. This difference could be 
a relic of the formation process with the current ortho-to-para 
ratio fixed at formation and/or set by the ratio of the parent 
species. Alternatively, it could indicate that not all of the 2i 2 - 
loi absorption seen at 107.7 /mi is associated with CH2, possi- 
bly due to blending with another line. Reducing the para col- 
umn density by approximately a factor of 2 would lead to an 
ortho-to-para ratio of 3. 

Overall, the column density summed over all velocities in 
both ortho and para states is N(CH 2 ) = (7.5 ± 1.1) x 10 14 cnr 2 , 
giving a final ratio with CH of [CH/CH 2 ] = 2.7+0.5. If the para 
column density is reduced to make the ortho-to-para ratio equal 
to 3, the [CH/CH 2 ] ratio would be increased to 3.7. 

5. Abundances 

To discuss CH and CH2 abundances, we need an estimate of 
the hydrogen column density towards Sgr B2. 

In the Galactic spiral arm clouds along the line of sight 
(-100 to +30 km s -1 ), HCN and CS measurements suggest 
average densities of 200 cm -3 , but also the presence of gas up 
to 10 4 cirr 3 (Greaves 1995). Their structure is probably similar 
to photodissociation regions (PDRs) with a molecular core and 
atomic skin at the surface, and UV illumination provided by the 
mean interstellar radiation field (Vastel et al. 2002). N(H2) has 
been estimated to be 5, 9, 4 and 14 x 10 21 crrr 2 respectively in 
the -100, -40, -25 and km s -1 features (Greaves & Nyman 
1996), with corresponding atomic hydrogen column densities 
approximately equal to 2, 4, 2 and 8 x 10 21 cm" 2 (Vastel et al. 
2002). In order to compare with diffuse cloud models, we cal- 
culate the abundances taking into account the total hydrogen 
particle column, N H = N(H) + 2N(H 2 ), although it is difficult 
to determine whether the observed CH and CH2 co-exists with 
both atomic and molecular hydrogen regions - they are proba- 
bly confined to one or more specific layers within each cloud. 



E. T. Polehampton et al.: FIR detection of methylene 



7 



Using our fitted column densities for CH in Table 3, we calcu- 
late an approximate abundance of N(CH)/N H ~(0.6-3)xl0~ 8 
for these clouds. According to our average [CH/CH 2 ] ratio, 
N(CH 2 )/N H is thus (0.2-1. l)xl(T 8 . 

The situation is more complicated at the velocity of Sgr B2 
as there may be contributions to the absorption from layers 
with widely different conditions. Comito et al. (2003) have 
found that a significant fraction of the water absorption to- 
wards Sgr B2 must be due to a hot (500-700 K), low den- 
sity layer (also observed in ammonia lines; Ceccarelli et al. 
2002; Hiittemeister et al. 1995) but that there must also be ab- 
sorption from the warm (40-80 K) envelope. Goicoechea & 
Cernicharo (2002) have derived a large OH abundance at the 
velocity of Sgr B2 suggesting that there is a strong UV field 
producing clumpy PDRs and a temperature gradient from 40- 
600 K through the envelope. The total H 2 column density in 
front of the FIR continuum has been estimated from the FIR 
spectrum to be (2-10) xlO 23 cnT 2 (Goicoechea & Cernicharo 
2001). Using this to calculate approximate abundances gives 
N(CH)/N H ~(0.5-2)xlO- 9 andjV(CH 2 )/jV H ~(0.2-0.7)xl0- 9 . 

6. CH 2 in other sources 

We have also searched the ISO Data Archive for observations 
towards other sources that may show the low-lying CH 2 lines. 

6.1. W49N 

Observations that cover the wavelengths of the CH 2 lines were 
carried out using the LWS L04 mode towards the active star 
forming region, W49 N. This is the strongest IR peak in the 
W49 A molecular cloud complex which is located at 11.4 kpc 
from the Sun and 8.1 kpc from the Galactic Centre (Gwinn 
et al. 1992). Its spectrum shows a strong thermal continuum 
in the FIR with a peak near 60 /im (Vastel et al. 2001). Line 
emission associated with the molecular cloud itself is centred 
at 8 km s _1 (e.g. Jaffe et al. 1987), but there are also several 
features due to intervening gas at velocities between 16 and 
75 km s- 1 (e.g. Nyman 1983). Vastel et al. (2000) have mapped 
the CO emission in the ISO LWS beam and found 7 velocity 
components in the range 35-70 km s _1 . Recently, Plume et al. 
(2004) have measured H 2 towards W49 and the same features 
show up in absorption. They arise because the line of sight in- 
tersects the Sagittarius spiral arm in two places (e.g. Greaves & 
Williams 1994). 

The strongest ortho-CH 2 ln-Ooo fine structure component 
at 127.6 fim (7=2-1) is clearly detected in absorption in the 
LWS FP spectrum. However, the spectral resolution and signal- 
to-noise ratio in the data do not allow a detailed fit using 
all the velocity components observed in CO and H 2 lines. 
Therefore, we performed a simplified fit using two components 
to represent the strongest features seen in the H 2 spectrum of 
Plume et al. (2004). We fixed the velocities (and FWHM) to be 
37 km s _1 (10 km s~') and 61 km s _1 (7 km s~') and found the 
best fitting optical depths after convolving to the LWS resolu- 
tion. Lower velocity components were not necessary to repro- 
duce the line shape, showing that all the absorption comes from 
the line of sight clouds. Only the 7=2-1 component was used 



: W49 N 
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Fig. 3. Data around the 7=2-1 (127.6 fim) and 7= 1-1 (127.8 pm) fine 
structure components of the ln-Ooo CH 2 transition toward W49 N. 
The solid line shows a model derived from a fit to the 7=2-1 compo- 
nent and the expected relative line strength of the 7=1-1 component. 

for the fit. The result is shown in Fig. 3, which also shows a pre- 
diction for the 7=1-1 (127.858 //m) fine structure component 
based its expected relative line strength. Column densities for 
the Ooo level, calculated from the best fitting optical depths and 
equation 2, are (0.8±0.4)x 10 14 ctrT 2 and (1 .2+0.4)x 10 14 cirT 2 
at 37 km s _1 and 61 km s _1 respectively. 

The para transition at 107 /zm (2 12 -l i) is not detected 
above the noise. Assuming that only the two lowest levels are 
populated and using a 2cr limit on the line depth, gives an ortho- 
to-para ratio >2.1. 

The CH line at 149 fim was not observed toward W49 N 
using the FP mode of the LWS, but it was measured with the 
lower resolution grating mode. At this resolution the A-doublet 
line splitting, as well as the line of sight structure, is not re- 
solved and we can only estimate a total CH column density 
that may include absorption associated with W49 itself, giv- 
ing 8.4 x 10 14 ctrT 2 . Sume & Irvine (1977) have observed CH 
radio A-doublet lines towards W49 A and they find column 
densities of 8 x 10 13 cirT 2 and 10 x 10 13 ctrT 2 in the veloc- 
ity intervals 30-45 km s _1 and 54-72 km s _1 . Rydbeck et al. 
(1976) also observed the A-doublet lines, giving column den- 
sities of 7 x 10 13 ctrT 2 and 18 x 10 13 cirT 2 . Both authors show 
that these two velocity intervals give the most important CH 
contributions. This means the extra column density seen in the 
LWS grating measurement probably also has highest contribu- 
tion from these ranges and the radio measurements may be an 
underestimate (probably due to the much larger beam size of 
15' used for these observations). 

Using the average of the two radio CH measurements and 
combining the ortho column density and para upper limit for 
CH 2 , gives a lower limit for the [CH/CH 2 ] ratio of 1.4 and 0.8 
in the two velocity intervals. 

6.2. Sgr A* 

Observations using the LWS FP mode were also made towards 
Sgr A* at the Galactic Centre. It has a FIR continuum show- 
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Fig. 4. The CH 2 Tli/2 7=3/2-1/2 transition observed towards Sgr A* 
showing two resolved A-doublet components. 



ing strong thermal emission from dust, peaking at ~50 yum. 
Observations of several hydrocarbons have been reported to- 
wards this source, including CH 3 (Feuchtgruber et al. 2000). 
We have searched LWS FP observations for the lines of CH2 
and CH and Fig. 4 shows the observed absorption due to the 
2 ni /2 7=3/2-1/2 (149 /an) rotational transition of CH. These 
lines show peak absorption near km s _1 , probably due to local 
gas. There is also absorption at negative velocities due to galac- 
tic spiral arm clouds. This picture is consistent with observa- 
tions of CH A-doublet lines towards Sgr A, which showed that 
the strong Galactic Centre component normally observed in 
other molecules at +50 km s _1 is relatively weak in CH (Genzel 
et al. 1979). At the wavelength of the ortho CH 2 1 11-O00 transi- 
tion (127 /j.m), the signal-to-noise ratio in the data is rather low 
compared to the Sgr B2 observations and we do not detect any 
absorption. The para 2i2-loi transition at 107.7 fim is also not 
detected. Combining 2cr upper limits for both transitions, we 
calculate a lower limit to the [CH/CH 2 ] ratio at km s _1 of 2.6. 

6.3. NGC7023 

The ISO database also contains LWS FP observations of CH2 
towards three positions in the NGC7023 PDR. However, no 
lines are detected above the noise in the spectra. This is con- 
sistent with LWS grating observations in which no OH, CH or 
CH2 were detected toward any position in NGC7023 (Fuente 
et al. 2000). 

7. Discussion 

7.1. Comparison with models 

In Table 4 we summarise the results from several chemical 
models for diffuse (van Dishoeck & Black 1986) and dense 
clouds (Lee et al. 1996), as well as a recent model that 
takes both gas-phase and grain-surface chemistry into account 
(Ruffle & Herbst 2001). In these models, the abundances of CH 
and CH2 are generally predicted to be higher towards lower 
densities and temperatures. We also include the results from a 



PDR model of the IC63 nebula from Jansen et al. (1995) - this 
is thought to have similar physical conditions to the Sgr B2 en- 
velope in ~ 10 5 crrT 3 ; 650-900 times the average interstellar 
UV field). In this model, CH and CH 2 peak towards the edge of 
the PDR at A v < 2, with [CH/CH 2 ] generally ~2 (although in 
a model with higher fractional abundance of carbon they cal- 
culated a value ~20). The CH abundance reached is relatively 
high compared to the dense cloud models with similar «h- In 
PDR models with higher density and more intense UV fields, 
even larger CH and CH2 abundances can be reached in the outer 
layers (Sternberg & Dalgarno 1995). 

Table 4 compares our estimated abundances with the mod- 
els. In the line of sight clouds towards Sgr B2 (-100 to 
+30 km s _1 ), Greaves & Nyman (1996) have examined the 
abundances of 11 species observed at 3 mm (HCO + , HCN, 
HNC, CN, CCH, C3H2, CS, SiO, N 2 H+, CH 3 OH and SO) and 
found that overall, the chemistry in these features is similar 
to the dark cloud TMC1. However, our abundances are much 
closer to diffuse cloud values, showing that CH and CH 2 may 
reside in the less dense parts of the clouds. In Sgr B2 itself, our 
derived abundances are much more uncertain due to the diffi- 
culty in estimating /V(H), but are lower than the line of sight 
clouds, consistent with denser material in Sgr B2. 

We can also compare our results with the previous obser- 
vations of CH2 in diffuse (Lyu et al. 2001) and dense clouds 
(Hollis et al. 1995). Using the results of Lyu et al. (2001) with 
hydrogen abundances for their lines of sight (see van Dishoeck 
& Black 1986; Snow et al. 1996), we calculate value consistent 
with our line of sight values and the diffuse cloud model listed 
in Table 4. For the high densities and temperatures prevailing 
in the Orion-KL and W51 M hot cores, Hollis et al. (1995) de- 
rive X(CH 2 ) ~ 4 x 10~ 9 and 7 x 10~ 10 , respectively, carefully 
considering beam-filling factors. These values are closer to our 
results for Sgr B2. 

7.2. CH/CH2 chemistry 

Both CH and CH2 are formed by the dissociative recombina- 
tion of CH3 and primarily destroyed by reaction with atomic 
oxygen to form HCO, a constituent of many more complex 
molecules. Hollis et al. (1995) compared their derived CH2 
abundance in the Orion-KL and W51 M hot cores with ra- 
dio observations of CH to determine branching ratios for the 
formation of CH and CH 2 of 0.1 and 0.9 respectively. The ac- 
tual branching ratios for the dissociative recombination of CHJ 
have since been measured in the lab by Vejby-Christensen et al. 
(1997). There are 4 important pathways, with CH2 as the major 
product (40%), two reactions producing CH (30%) and the re- 
mainder producing atomic carbon (30%). This appears to agree 
with models of diffuse and translucent clouds, which required 
a high abundance of CH2 to reproduce the measured column 
densities of C 2 (van Dishoeck & Black 1986, 1989) 

The above discussion indicates that CH2 should be more 
abundant than CH. However, our observations show that in 
the line of sight towards Sgr B2, CH2 has less than half the 
abundance of CH - if this were not the case, we would have 
detected much stronger CH2 absorption. The observations to- 
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Table 4. Comparison of model predictions for diffuse and dense clouds with the observed abundances (compared to 2iV(H 2 )+N(H)) of CH and 
CH 2 (X(CH) and X(CH 2 )). The diffuse cloud model is for the line of sight to f Per (van Dishoeck & Black 1986): vDB86. The quoted density 
and temperature are for the centre of the cloud in the model. The dense cloud models are steady state values corresponding to the gas phase 
"new standard model" of Lee et al. (1996): L96 and the gas-grain model of Ruffle & Herbst (2001): RH01 including photochemistry (Pl/A at 
10 8 yr). The PDR model is for a density of 10 5 cirT 3 and radiation field -650-900 times the average interstellar value (Jansen et al. 1995): J95. 
The observational results derived from our ISO measurements are shown as well as previous results from Lyu et al. (2001) and Hollis et al. 
(1995). 
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0.039 


0.16 




Dense (RH01) 


2x 10 4 


10 


0.1 


1.2 


0.08 




PDR (J95) 


10 5 




1.4 


0.7 


2.0 


Results from 


Sgr B2 Line of Sight 






0.6-3 


0.2-1.1 


2.7±0.5 


this work: 


Sgr B2 
W49N37kms-' 
W49N61 kms-' 

SgrA' 






0.05-0.2 


0.02-0.07 
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£Oph 
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4.3" 
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results: 
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100-200 


0.06 c 
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0.1 f 




W51 M 


~ 10 4 


100-200 


0.007 f 


0.07" 


0A e 



" Combining A'(CH 2 ) from Lyu et al. (2001) with a hydrogen column density as quoted in van Dishoeck & Black (1986). 
* Based on N(CH)=2.5xl0 13 cnr 2 from van Dishoeck & Black (1989). 

c Combining A^CIL) from Lyu et al. (2001) with a hydrogen column density from Snow et al. (1996). 
d Based on N(CH)=8xl0 13 cnr 2 from van Dishoeck & Black (1989). 

e Taking values as quoted by Hollis et al. (1995) - but note that X(CH) could be significantly higher - see text. 



wards W49 and the non-detection towards Sgr A* also appear 
to confirm this. Furthermore, the results of Lyu et al. (2001) 
towards HD 154368 combined with previous measurements of 
CH (see Table 4) show very good agreement with the Sgr B2 
ratio. 

These results are contrary to the previous observations to- 
wards Orion-KL and W51 M where CH2 was found to be 10 
times more abundant than CH. This may be due to the difficulty 
in determining the CH column density equivalent to the mea- 
sured CH2 emission in these sources. Hollis et al. (1995) used 
CH column densities derived by Rydbeck et al. (1976), giving 
N(CH) = 2 x 10 13 cirT 2 towards W51 M. In contrast, Turner 
(1988) made a detailed study of the A-doublet satellite lines 
from the first rotationally excited state towards W51, giving a 
value for the 57 km s _1 velocity component that is more than 
100 times larger: 6.2xl0 15 cm 2 . This indicates that [CH/CH 2 ] 
could be as high as 30, although the radio CH line observations 
are likely to sample a much larger and cooler region than the 70 
GHz CH2 lines which require high temperatures to be excited. 

This shows it is extremely important to have consistent 
measurements of CH and CH2 (with transitions at similar ener- 
gies, using the same beam size and with consistent calibration) 
in order to be able to make a good comparison. So far the cur- 
rent data towards Sgr B2 provides the best comparison of the 
two species. 



In order to explain the low abundance of CH2 with re- 
spect to CH observed towards Sgr B2, we require other for- 
mation/destruction processes to be important for the CH/CH 2 
balance. The PDR model of Jansen et al. (1995) shows that the 
presence of strong UV radiation can reproduce our observed 
[CH/CH2] ratio. This may explain the ratio in Sgr B2 but in 
the galactic spiral arm clouds, the UV field is much weaker 
and densities lower than used by Jansen et al. (1995). Viti et al. 
(2000) have modelled the hydrocarbon chemistry in diffuse and 
translucent clouds (n=300 cirT 3 ) including the interaction of 
gaseous C + with dust grains. This leads to the production of hy- 
drocarbons by surface reactions. In some of their models (for 
A v =4), they find column densities approaching those that we 
find in the line of sight clouds towards Sgr B2 and [CH/CH2] 
ratios 7.7-13.6. These surface reactions could explain our high 
observed [CH/CH 2 ] ratio. 

8. Summary 

We have made the first detection of the low-lying rotational 
transitions of the CH2 molecule in the ISM towards Sgr B2 and 
W49 N in absorption. We do not detect CH 2 towards Sgr A* or 
NGC7023. For Sgr B2 we were able to compare these obser- 
vations with measurements of the ground state rotational lines 
of CH, providing a good estimate of the total column density 
of both species along the line of sight. These observations pro- 
vide the best comparison of the two species to date, giving a 



10 



E. T. Polehampton et al.: FIR detection of methylene 



[CH/CH 2 ] ratio of 2.7+0.5, probably fairly constant in all the 
observed velocity components. The results towards W49 N ap- 
pear to agree with this ratio, giving a lower limit of ~ 1 .0. 

Comparison with chemical models shows that the abun- 
dances in the line of sight clouds are close to diffuse cloud 
conditions whereas in Sgr B2 itself they indicate denser gas is 
present. Our high [CH/CH2] ratio can be explained by models 
including grain surface reactions (e.g. Viti et al. 2000). 

Future observations in the FIR with SOFIA and in the 
sub-millimetre with telescopes such as the Atacama Pathfinder 
Experiment (APEX) will be highly interesting to extend the 
study of CH2 to other sources. 
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Appendix A: Observations used 



Table A.l. Log of the observations used. 
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* Blended with NH line. 



